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Abstract 



^ ■ A study of prompt photon plus associated jet production is performed at next- 

to-leading order (0(aa 2 )) in QCD at ^/S = 200 — 500 GeV, appropriate for 



the RHIC polarized collider experiment. Momentum correlations between 



ON I the jet and photon are examined and the utility of the process as a method for 

constraining the size and shape of the polarized gluon density of the proton 
Q\ • AG is examined. 



ST! I- INTRODUCTION 



^. . 

With the advent of RHIC, the Relativistic Heavy Ion Collider at Brookhaven, QCD will 
^ ■ enter into a new interesting phase in which polarized high energy collisions will become a 
standard tool of analysis in high energy physics. RHIC, as a pp collider, will be endowed 
with a very high luminosity (500 pb -1 ), not easily accessible to pp colliders and will be 
spanning a center of mass energy range between 50 to 500 GeV. 

One of the main programs at RHIC will be to nail down the size and the shape of the 
polarized parton distributions which, at the moment, suffer from significant model depen- 
dence, especially in the gluon contribution (AG). Interest in spin physics has grown in the 
last few years thanks to the various polarized DIS experiments and to the existing proposal 
for the construction of a pp collider at DESY (HERA-iV) running at lower energy (^ 50 
GeV). 
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More generally, besides the important information on the spin structure of the nucleon 
which can be gathered from these experiments, it is obvious that one might try to look 
for physics beyond the Standard Model using polarization as a tool to suppress unwanted 
background and to enhance specific signals. The studies of chiral couplings in the Standard 
Model or even compositeness will require precise measurements of the hadronic background. 
There have been many attempts in the last few years at creating a backbone-program for 
RHIC by analyzing a set of processes in leading order and by defining suitable observables 
which are more easily accessible to the experimental investigations. 

More recently, next to leading order (NLO) studies of various processes (see for instance 
]1|]) have been presented. A partial review of these developments can be found in 0. Among 
the most interesting processes which require an accurate theoretical determination are single 
and double prompt photon production, single and 2-jet production and the Drell-Yan lepton- 
pair production. Work on the Drell Yan cross section at NLO has been presented 
limited at the moment to the non-singlet sector, which is sensitive to the polarization of 
the quark distributions. Due to the fact that the initial partons are longitudinally (or 
transversely) polarized, the calculation of the hard scatterings are far more involved than in 
the unpolarized case. 

Although the studies of the total cross sections for some of these and other related 
processes will be necessary in order to guarantee experimental observability of the cross 
sections and thus the spin asymmetries, which are usually estimated to be small, the study 
of the event-structure of these processes can provide interesting new insights into the spin 
distributions. In this paper we study and compare the structure of prompt photon plus single 
jet production in the polarized and unpolarized cases at RHIC center-of-mass energies. 

II. MONTE CARLO CALCULATIONS AT NLO 

Compared to total cross sections, distributions usually reveal more detail about the 
underlying hard scattering mechanisms as well as more information on the x-dependence 
of the model structure functions. The drawback is that they are generally more difficult 
to define theoretically to NLO, since many of them are affected by non-canceling infrared 
divergences. Modern developments in combined Analytic/Monte Carlo techniques allow to 
get exact numerical results for the NLO corrections in a reasonable amount of time. In 
addition they allow great flexibility in placing experimental selections such as jet definitions 
and photon isolation cuts on the cross sections, thereby allowing a more realistic and direct 
comparison with data. 

The analytical part of the calculation involves the 1) exact (analytical) evaluation of 
the virtual corrections and 2) the exposure (by a cutoff regularization) of all the mass 
singularities in the real emissions. The rest of the phase space is then integrated over 
numerically. We omit a general presentation of the method which can be found elsewhere 
and focus on the study at NLO of the rapidity correlations in photon plus associated jet 
production. 

In the next sections, after a brief overview of the various contributions to the cross 
section, we move to a study of the rapidity correlations between the photon and the jet. 
The analysis presented is similar in spirit to the study of momentum correlation given in ref 
. with the due modifications. 
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III. THE PHOTON PLUS JET CROSS SECTION 



The main goal at polarized hadron colliders is the study of the polarized parton distri- 
butions of the nucleon which are defined by 

A/ 4 (x, M 2 ) = f+(x, M 2 ) - fr(x, M 2 ). (3.1) 

The corresponding unpolarized distributions are defined by 

f^M 2 ) = f?(x,M 2 ) + f-{x,M 2 ), (3.2) 

where and f~ represent the distribution of partons of type i with positive and negative 
helicities respectively, with respect to that of the parent hadron. The hard subprocess 
scattering cross sections are defined by 

A<7y = l - (a ++ - o + -) (3.3) 

and 

&* = \ (^ ++ + ° + -) ■ (3-4) 

for the polarized and unpolarized cases respectively. 

It has been observed by many authors that the cross section for prompt photon produc- 
tion is dominated by the subprocess qg — > ^q in hadronic collisions already in leading order. 
This means that the cross section, if properly understood, could potentially prove very useful 
for providing information on the unpolarized gluon densities, g(x, Q 2 ), of hadrons. Original 
Born level studies of this cross section also indicated that the same is true in the polarized 
case. It has therefore been suggested that it may prove useful in pinning down the polarized 
gluon densities 0. In this context it has been examined quite extensively in leading, and 
more recently in next-to-leading order [g] [|]. The most recent NLO study |L(J included the 



effects of photon isolation on the cross section. The results of all the NLO studies confirmed 
the conclusions from the LO ones that the cross section is very sensitive to AG and that 
the asymmetries are perturbatively stable. 

Recently the photon plus jet cross section was studied in LO [O] and NLO [T^| and 



HERA-N cms energies. The main conclusions from these studies is that this cross section 
will give more detail about the x-shape of the polarized gluon distribution than the single 
prompt photon cross section. The present calculation follows along the same lines as that in 
ref. [0] and we therefore do not give more detail about the contributing subprocesses here, 



except to say that the fragmentation contributions are estimated in LO here as well. In the 
present case, before isolation, these contributions are numerically more important due to the 
higher cms energy of the RHIC collider. The results of ref. [IIJ also showed that isolation 



significantly reduced these contributions and that that their presence, although affecting the 
predictions for the cross sections, did not affect the asymmetries very much. 
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A. Rapidity Correlations to LO and NLO 



As mentioned above, although single inclusive prompt photon production will definitely 
be very important for constraining the size of AG, information on the detailed rr-shape of the 
distribution will not be as easily extracted. This is because the calculation of the inclusive 
cross section involves one convolution over the momentum fractions, x, of the initial partons. 
In fact, at hadron level, the factorized hadronic total cross section is generically denoted as 

&o = Y,f dx 1 dx 2 Af i (x 1 , )Af j (x 2 ) J dAa. (3.5) 

where we sum over all the partons i,j. The practical effect of the integral over the Xi is 
that a measurement of the kinematic variables of the photon is not sufficient to determine 
them. If, on the other hand, one or more of the jets produced in the reaction is also tagged, 
no convolution is involved in the calculation and the cross section is directly proportional 
to the parton densities. 

The cross section of interest here is the triple differential cross section 

d 3 Aa^ J 



dpTpdrfi dr\ 3 ' 



where rp and r/ J are the pseudorapidities of the photon and jet respectively and is the 
transverse momentum of the photon. 
We use light-cone coordinates 



Pi = p{n = xiQn p 2 =p 2 n = x 2 Qn n = l/v2(l, _L, ±1) 

with Q = \/~S /2 denoting the ± components on the light cone of the two incoming hadrons 
of momenta P\ and P 2 . We also set p\ = x\Pi and p 2 = x 2 P 2 for the two partons that enter 
the hard scattering. 

The 4-dimensional 5 function and the integration variables ki are also rewritten on the 
light cone (ki = kf , k^ki^) and after simple manipulations we get 

Aa ~ E f dx ^ f dxzAfiMAfjfa) J dktdk 2 d 2 k l± 5{kl)5{kl)A\M\ 2 . (3.6) 
We have set 



kf = ^ k 2 = ^e--. (3.7) 



We have defined the two rapidities 



y, = ^log^ < = 1,2 (3.8) 



and introduced the rapidity difference Ay — yi — y 2 - 
We easily get 
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d 3 Aa^ J dAa , ^ A „, SA , , , dAa , , 

2A; ± xTxi^A/ i (3Jr)A/ i (^)— - (3.9) 



d'Pj.dvpdr] 3 dk±dyidi/2 ^ dt 

A derivation of this result is illustrated in the appendix. A similar result applies to the 
unpolarized case, except that the polarized cross sections and structure functions are re- 
placed by the corresponding unpolarized ones. In this next-to-leading order calculation a 
jet definition is required. Throughout we use the Snowmass |18[ jet definition. 

All the 2-to-2 contributions (Born and virtual) to the photon-plus-jet cross section give 
contributions with structure functions sampled at fixed kinematical points x,. Thus the 
double longitudinal asymmetry A LL , defined as the ratio of the polarized to the unpolarized 
cross section, is directly proportional to the ratio AG/g JTTJ] in kinematic regions where 
other subprocesses such as qq scattering can be neglected. This guarantees sensitivity of the 
asymmetries to AG. 



IV. RESULTS 



All results are displayed for pp collisions at the center-of-mass energy = 200 and 500 
GeV which are energies typical for the RHIC experiment at Brookhaven. For the unpolarized 
cross section the CTEQ4M parton densities |H| are used throughout, and the value of A^j^ 



corresponding to this distribution is also used. Use of other unpolarized parton densities at 
the x-values probed here do not yield significantly different results. For the polarized case 
the GRSV and GS (15| distributions are used with the corresponding values for Ajjg. 



The authors of ref. [14] and [[HJ have proposed various parameterizations of the polarized 
parton densities differing mainly in the choice of input for the polarized gluon density AG. 
In the case of the GRSV distributions we use the 'valence' set which corresponds to a fit 
of the available DIS data (referred to by the authors as the 'fitted' AG scenario), the large 
gluon fit which assumes that AG(x, Ql) = g(x, Ql) at input (the 'AG = g' scenario) and the 
small gluon fit which uses AG(x, Qq) = at the input scale (the 'AG = 0' scenario), which 
in this case starts at the very low value of Qq = 0.34 GeV 2 . The latter two distributions are 
intended to represent extreme choices for AG. These parameterizations give gluon densities 
which differ in their absolute sizes as well as in their x-shape. The GS parameterizations 
provide three fits to the data; GS A, GS B and GS C. It has been shown that the GS A and 
GS B distributions do not differ very much from the AG = g and fitted AG sets of GRSV 
respectively, whereas the the GS C set is widely different from any of the others. We shall 
present distributions using the three GRSV sets discussed above, along with the GS C set for 
comparison. For the fragmentation functions we use the LO asymptotic parameterization of 
ref. ||16|| . As will be shown, the choice of fragmentation functions makes very little difference 



to the predictions, since these processes account for only a small fraction of the cross section. 

The renormalization, factorization, and fragmentation scales are set to a common value 
fi = Pt unless otherwise stated. Since there are two 'particles' in the final state, the jet 
and the photon, both of whose transverse momenta are large, an alternative choice might 
be fi = p 3 T or some function of p\ and p 3 T . The results of the calculations show that the 
magnitudes of pTp and pj> tend to be comparable and that dependence of the asymmetries on 
fi is slight, although the individual cross sections may vary significantly with \i. Therefore, 
choices of \x different from fi = p\ should not produce significantly different predictions 
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for the asymmetries . The two loop expression for ot s {Q 2 ) is used throughout, with the 
number of flavors fixed at Nf = 4, although the contribution from charm was verified to 
be negligible at the energies considered and is not included. None of the currently available 
polarized distributions include a parameterization of the charm quark distribution. A new 
NLO parameterization is in preparation which includes charm, but is not yet available for 
this study Finally, the values of both the jet cone size and photon isolation cones are 
fixed at Rj = 0.7 and r 7 = 0.7 respectively, unless otherwise stated. 

Figs. la and b show the triple differential cross section as a function of p\ of the photon 
for the various parameterizations at v^S* = 500 and 200 GeV respectively. The unpolarized 
cross sections are shown for comparison. The curves were obtained by averaging over bins 
Apj> = 1 GeV and the restriction pj* > 10 GeV was imposed. In addition both the pho- 
ton and jet rapidities are averaged over the central region, —0.5 < r]' y ,r] J < 0.5. All the 
parameterizations give distributions which are distinctly different in both their shapes and 
sizes. Their relative sizes are in direct relation to the sizes and shapes of their respective 
gluon distributions. This is most obvious for the GS C parameterization which has the most 
distinct gluon distribution, being negative over part of the x-range. The curves show a rise 
between p\ = 10 and 11 GeV because of the restriction p^, > 10 GeV and the fact that the 
bin centered at 10 GeV is averaged over the range 9.5 < p^ T > 10.5 GeV. Above p\ = 10 GeV 
both the two- and three-body contributions to the cross section are finite whereas below this 
value, only the latter are finite since the two-body contributions always produce a photon 
balancing the px of the jet. 

Figs.lc and Id show the asymmetries for la and lb respectively. There are clear distinc- 
tions between the predictions for the various parameterizations which will certainly make 
them distinguishable in the experiments. As expected, at the lower the asymmetries are 
larger, but the corresponding cross sections are smaller. Combining both sets of results will 
nevertheless suggests that a measurement of the polarized gluon distribution will be possible 
between x = 0.04 and x = 0.5. 

In Figs. 2a and 2b we look at distributions in rp at p\ = 10 GeV with the restriction 
p J T > 10 GeV still imposed. In both cases yS = 200 GeV. In Fig. 2a the jet is restricted to be 
in the central rapidity region, —0.5 < r] J < 0.5, and in Fig.2c it is restricted to the forward 
region, 0.5 < r] J < 1.5 . The textures of the different curves follow the same conventions 
as in figs.la-d. A visual comparison of the figures show clearly that when the jet is in the 
central rapidity region the, the photon rapidity distribution peaks at 77 7 = 0, whereas when 
the jet is restricted to the forward rapidity region, the 77 7 also peaks in this region. This 
positive rapidity correlation between the photon and jet is present for both the polarized 
and unpolarized cross sections although it is clearly stronger in the unpolarized case. 

Fig. 2c which shows the corresponding asymmetries for the curves in fig.2a, verifies that 
the distributions are symmetric and also shows that they have distinctive shapes. The 
corresponding asymmetry curves for fig. 2b all rise sharply in the negative rapidity region, 
i.e., in the opposite direction to the region where the rapidity distributions peak. Thus 
the asymmetries display a negative rapidity correlation between the photon and jet. An 
explanation for this behavior was given in ref. ]12|] in terms of competing effects between 
the subprocess matrix elements, which tend generate positive rapidity correlations, and the 
polarized parton distributions which tend to produce negative correlations. 

In figs. 3a and 3b rapidity distributions are plotted for the jet using similar cuts to those 
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used in figs. 2a and 2b. In fig.3a, the photon is restricted to be in the central rapidity region 
whereas in 3b it is restricted to be in the forward region. As before, there is a net positive 
rapidity correlation between the photon and jet for the polarized as well as the unpolarized 
cases, but as fig.3d shows, the effect is stronger for the unpolarized case. This leads the 
asymmetries to peak at negative rapidities of the jet. 

One of the striking features of the asymmetry curves in figs. 3c and 3d is the differences 
in their shapes as compared to each other and as compared to those of figs. 2c and 2d. The 
parameterizations with the larger polarized gluon distributions the "fitted" AG and AG = g 
scenarios give asymmetries which decrease as r\ J moves away from the central region in fig. 3c 
and away from r\ J ~ —1.5 in fig.3d. The differences between these curves and those of fig. 2 
are explained by the asymmetric pr cuts between the photon and jet which affects the phase 
space available for jet or prompt photon production differently. 

The differences between the shapes of the asymmetries for the various models of polarized 
parton distributions in figs. 3c and 3d are due to the differences in the x-shapes of the 
polarized parton distributions, particularly those for AG. All the asymmetries fall as i] J 
goes to the extreme values as may be expected since the unpolarized gluon distribution 
has a singular (x~ a where a < 1) behaviour at small x, whereas all the polarized gluon 
distributions go to zero at small x. All the other differences between the asymmetries are 
due directly to the differences between the corresponding polarized parton distributions. 
This re-enforces the conclusion that this cross section will undoubtedly yield very important 
information on the polarized parton distributions. 

V. CONCLUSIONS 

We have examined the possibility that both the size and x-shape of the polarized gluon 
distribution of the proton, AG, may be measured at RHIC via a measurement of the photon 
plus jet cross section. Control over the kinematic variables of both the photon and jet allows 
a much better determination of the rr-value probed when compared to inclusive prompt 
photon production. A comparison of the predictions obtained using different polarized 
parton densities show that a clear distinction between both the sizes and shapes should 
be possible. Further more detailed information on the x-shapes of the polarized gluon 
distribution will be obtainable by carefully choosing the kinematic regions in which the jet 
and/or photon is tagged. 

Assuming that the 'fitted AG' scenario is the most plausible distribution, then a typical 
value for the asymmetry, All is 5%, but given the uncertainty in AG the asymmetry could 
be as small as 1% or as large as 20%. The expected small- a; behavior of the polarized 
and unpolarized distributions lead to differencs in correlations between the rapidities of the 
photon and jet in each case. This effect can produce very large asymmetries in certain 
regions of phase space which may be expolited and used to discriminate between various 
models of AG. 

VI. APPENDIX 

The cross section for the 2-to-2 process is generically given by 
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\M\ 2 (2tt)\ 4/ , , 5 + (k 2 )\ 4/ 

= ' ^ S 4 (p 1 +p 2 -h- k 2 ) +l y +) 2 ^ 5 A ( Pl +p 2 -h- k 2 ), 6.1 



and working in the cm frame 



' |M| 2 . (6.2) 



dt 16vr 2 

We expand the final state momenta in the light cone variables 

cr ~ J dkfdk^d 2 k 1 j_dk2dk 2 ~d 2 k 2 ±5 ( - 2 \ki± + k 2 ±)5(pf + p 2 — k+ — k£) 
x5(pi+p 2 - k± - k 2 )5+(k 2 1 ) 5+(k 2 2 ) 

and integrate over k±, k 2 and k 2 ±, thereby eliminatigng 3 of the 5 delta functions. 
The remaining 2 delta functions are rewritten as 

.2 ' 



5{k\) = 5 (V2k ± e Vl (x 2 Q - k ± e~ y2 ) - k] 
5{k 2 2 ) = 5 (v^e-^iQ - k ± e yi ) - k] 



and the integration over the parton fractions x±, x 2 performed by the relation 



with 



x 2 



y/2Qe* 
k ± {l + e A y) 
V2Qe-y* 



Therefore, the convolution integral for the cross section is reorganized as follows 

da~ ( dxi f 1 dx 2 f(x 1 )f(x 2 ) [ dkfdk 2 d 2 k ± 5(k 2 1 )5(kl)\M\ 2 
Jo Jo J 



(6.3) 



(6.4) 



5+(kl) 5 + (k 2 2 ) = 2k 2j A yQ2 6 ( Xl - X l) 5 ^ - X 2) ( 6 ' 5 ) 



k±(l + e A y) 

Xi = 



(6.6) 



dxi dx 2 f(x 1 )f(x 2 ) J dkjdk 2 d k ± 2fc2 6 a v q2 \ M \' 

At this point we change the remaining integration variables to the rapidity space using 
d(kf, k 2 )/d(yi,y 2 ) = exp[Ay]k 2 L /2 and perform the integration over one of the angle in k±. 
With the correct normalization, we get the result in |3.9|. 
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FIGURE CAPTIONS 



[1] (a) p\ distribution of the photon plus jet triple differential cross section 
d 3 a" /J /dp' y T dri"'dri J at y/S = 500 GeV for various polarized parton distributions 
and for rapidities of the photon and jet averaged over the region —0.5 < rf' , i] J < 
0.5. The cut p J T > 10 GeV is imposed, (b) Similar to (a) at VS = 200 GeV. (c) 
and (d) Corresponding asymmetries for the distributions in (a) and (b) respec- 
tively. 

[2] (a) Distribution in the rapidity of the photon at p? = 10 GeV and i] J averaged 
over the region -0.5 < 7] J < 0.5 and p J T > 10 GeV at = 200GeV. (b) Similar 
to (a) but for 0.5 < i] J < 1.5. (c) and (d) Asymmetries for the curves shown in 
(a) and (b) respectively. 

[3] (a) Distribution in the rapidity of the jet at p\ = 10 GeV and rp averaged over 
the region -0.5 < ^ < 0.5 and p-{ > 10 GeV at v 7 ^ = 200GeV. (b) Similar to 
(a) but for 0.5 < rp < 1.5. (c) and (d) Asymmetries for the curves shown in (a) 
and (b) respectively. 
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